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We suggest a systematic method of extension of quasi-exactly solvable (QES) 
systems. We construct finite-dimensional subspaces on the basis of special functions 
(hypergeometric, Airy, Bessel ones) invariant with respect to the action of differential 
operators of the second order with polynomial coefficients. As a example of physical 
applications, we show that the known two-photon Rabi Hamiltonian becomes quasi- 
exactly solvable at certain values of parameters when it can be expressed in terms 
of corresponding QES operators related to the hypergeometric function. 
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I. INTRODUCTION 

Usually, the models of the quantum mechanics are considered to be exactly-solvable (when 
the eigenvalues and eigenfunctions are known) or non-solvabe at all (when the eigenvalues 
and eidenfunctions are unknown, and they can be found numerically only). It became a 
great surprise that an intermediate case is also possible, when inside the Hilbert space there 
is an invariant subspace for which the eigenvalues and eigenfunctions could be found from 
algebraic equations. This type of systems is called quasi-exactly solvable (QES). For the 
wide range of the one-dimensional QES the Hamiltonian possesses hidden Lie algebra sl(2) 
and represents a bilinear form of first-order differential operators which satisfy the same 
commutator relations as the spin ones [1], [2], [3], [4]. The QES Schrodinger operators based 
on the s/(2) representation were studied in [1], [5], [6], [7], [8]. Later it became possible to go 
beyond the Lie algebraic context. In particular, starting from the second-order differential 
operators 
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such that they possess the finite-dimensional representations both in the subspace V n = 
span{l, x, a;™ -1 , x n } of all monomials of the degree < n and in the monomial subspace 
span{l, x 2 ,... , x n ~ 1 ) x n } [9], [10], [11]. In [12], [13], .the extension of the the corresponding 
subspace was suggested 

v = V n U g(x)V n = span{l, x 2 , x n ~\ x n , g(x),g(x) ■ x, g(x) ■ x n } (2) 

for certain functions g(x) [12], [13]. 

In the present paper we consider the problem of constructing QES operators which pre- 
serve subspaces of a more general form: 

M n = span{f\(x), f 2 (x), fn-i(x), fn(x)}. (3) 

The general features of differential operators with invariant subspace M n (3) were discussed 
in [15]. The main result obtained is the following: any linear operator P n G V(M n ), that 
preserve M n (3), may be represented as follows (a formula 5.8 in [15]): 

where are arbitrary constants, R u G V are arbitrary linear differential operators, V = 
V (J 7 ) is a space of the differential operators with coefficients belonging to the functional 
space T (hereafter denoted Vn)-, K m are operators belonging to the annihilator A (M n ): 

A (M n ) = {KeV\K (f(x)) = 0, for all f(x) G M n }, (5) 

Lj— are operators belonging to the affine annihilator /C (M n ): 

K (M n ) = {LeV\L (f(x)) = c G R, for all f(x) G M n } (6) 

c is a constant. All the definitions and symbols are in agreement with the ones used in [15]. 

In spite of the fact that the theorem of existence specifies a general form of the required 
differential operators (4), their finding is, generally speaking, not trivial task. The purpose 
of the given work consists in explicit construction of such QES operators, for which the set 
of functions f n (x) does not reduce to polynomials and represents special functions. With 
this purpose, we offer a method of QES extension. It consists in the following. To construct 
the QES Hamiltonians on the basis of the subspace (3), we select in (6) the subspace of 
differential operators K, 2 (M n ) of the degree two or less: 

/C 2 (M n ) = {L G K (M n ) I order(L) < 2}. (7) 
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The general approach suggested in [15] does not supply us with information, whether or not 
the space /C 2 {M n ) (7) is empty, without the constructing K, {M n ) for the chosen M n (3). 
To simplify matter, we start from the two-dimensional subspace M2 = span{fi(x),f2(x)}, 
for which operators Li e /C 2 {M2) are known to exist and can be explicitly constructed [15]. 
Let us try to extend the subspace M. 2 by adding function fz{x) to it in such a way that 
the new subspace M3 = span{fi{x), f 2 {x), f 3 {x)} be the three-dimensional and the space 
of corresponding operators /C 2 {M3) (7) be non-empty. Thus, transformation M 2 — ► M3 
has to lead to the change of the coefficients at in the operators L { e /C 2 {M.2), with the 
order of the operators e /C 2 (A / ( 2 ) preserved: order {Li) < 2. Repeating this operation 
n times, we obtain the QES-extension of .M 2 . In other words, the initial subspace M2 = 
span{f\{x), / 2 (a;)} permits the QES-extension if M2 may be extended to M n under the 
condition that the affine annihilator /C 2 {M n ) (7) is not empty. (It is worth noting another 
method of constructing invariant subspaces (3) which is based on the conditional symmetries 
[14], but it will not be discussed in the present paper.) 

II. CONSTRUCTION OF THE INVARIANT SUBSPACES 

Below we suggest the method that enables one to construct invariant subspaces M n . 
Let us start from the simplest case n — 2. Let us choose a linear independent basis 
{fo~{ x )i fo{ x )} f° r t ne subspace M2 = span{/ l "(a;), fo{x)} in the following way. Let 
us consider the function f{x), which satisfies the second order homogeneous differential 
equation with polynomial coefficients: q{x)-^f{x) + p{x)-^f{x) +r{x)f{x) = 0. The case 
when f{x) = constant ■ rr:f(x) is eliminated. We select the function f${x) = f{x) and its 
derivative f${x) = ^ f{x) = f'{ x ) as the basis components of 7W 2 : 

M 2 = span{f£(x), f {x)}. (8) 

For simplicity, we shall work directly with the operator (4) instead of the operators Li 
representing affine annihilator /C {M2) (6). 

Our strategy can be described as follows: 

I We find a general form of the operator of the second order P 2 for which subspace M2 
(9) is preserved. The existence of such an operator is supported of Lemma 4.10 [15]. 
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II We make extension of the subspace M 2 — ► M-a = span{ft, ft , fo , fi }• 

III We find a general form of the operator of the second order P4 for which subspace M.4 
is preserved. If a non-trivial solution P 4 7^ const exists for the chosen way of the 
extension, we pass to item IV, otherwise we change a way of extension. 

IV We make comparison of two operators P 2 and P4 and try to guess a general form of 
coefficients which enable us to repeat extension in the chain ( M 2 — > M.4 — > ... — > 
•M2(iv+i)) (V = 0,1,2,...) ) that leaves these subspaces invariants. As a result, we 
obtain the explicit form of operator P2(v+i) that acts on the elements of the subspace 
M 2 ( N +i) = {ft, ft, ft- ft, fo, fi, f 2 , ■-, }, N = 0, 1, 2, .... This is not the end 
of story since the general form does not fix by itself the concrete expression for the 
coefficients. 

V With the guessed general form of the operator P 2 (n+i) at hand, we require that it 
leave the corresponding subspace invariant, whence we find the concrete values of its 
coefficients. 

For explicit demonstration of the aforementioned algorithm, we shall consider a concrete 
choice of the function ft(x) : ft(x) = Fi [~',x] [16] where F 1 [~',x] is a hypergeometric 
function, fo(x) = F 1 [ s ~ 1 ;x] is its derivative up to a constant factor. 

The differential equation which the hypergeometric function obeys 
( x ' 3^ + s ■ — ij fo(x) = as well as its other properties can be found, e.g., in 
[16]. Then, we have: 

M 2 = span { F, [-;x], F 1 ;x]}, (9) 

Let us choose the operator of the second order P 2 = a 3 (x)-^ + o, 2 (x)-^ + a (x) and write 
down the condition of invariance of the subspace M. 2 (9) for the given operator: 

P2 (/o + ) = ci • fo + + c 2 ■ fo (10) 

P2 {fo) = c 3 • fo + C 4 • fo (11) 

Using rules of differentiation £ f + = & - £ f = s (s + 1) ^ + (s + s 2 + x) ^ and 
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equating coefficients at functions fg, f (10, 11) we obtain the solution: 



a 3 (x) = -c 3 • — + c 2 • s ■ x, 

s 

a 2 (x) = -c 3 • x + c 2 • s (s + 1) , (12) 

Oi (x) = c 3 • - - c 2 ■ s + c 4 + c 3 

s 

Ci = c 4 + c 3 . Substituting obtained coefficients a k (x), (k = 1,2,3) (12) into the expression 
for the operator P 2 we have: 

P 2 = c 2 -s-P 2 l ---P 2 2 + c 4 + c 3 + s-c 2 (13) 

s 

where 

tf = *£ + (« + i)-£.tf = ^ + «-4-*- (14) 

One can check that operators P2 1 , P 2 2 act on the basis functions according to the formulas 

P\ Uo) =fo+ l s - fo, Pl Uo) = fo, Pi (/o + ) = 0, Pi {Jo) =s-f -s-f + and, thus, 
preserve the subspace Ai 2 (9). This is agreement with the fact that the affine annihilator 
K, 2 (M) has the property of the vector space and possesses the operator basis (14) (see in 
[15]). 

The extension of the subspace M 2 (9) can be realized in many ways. The basic require- 
ment here consists in that after adding new elements of the basis M. 2 U span {/{*", f{} = 
span {fo, fi, fo, fi} = Ma the operator basis of P 2 (13, 14) should remain two-dimensional 
for the new subspace M.^. We consider the simplest way of extension - multiplication of the 
basis functions by the power function x n : M. 2 — > M.^ M. 2 (n+i) {N = 0, 1, 2, ....), 

{/o") fo } ~ > {fo ' /i + > fo -i fi } — *■ ••• ~~ > {foi A + ' •••/iV' fo ' /i ' /at} (15) 

where 

fn =X ' fn-1 = X ' /o") fn =X ' fn-1 = X ' fo ) (^) 

n = 0, 1, .., N — 1, N. To demonstrate that suggested extension is suitable, we shall consider 
subspace .M4 = {/q 1- , /q", A~} and shall write down the condition of its invariance with 
respect to the operator P 4 = b 3 {x)^ + b 2 (x)-^ + b (x): 

P±(j)=d-J (17) 

where / = {/ + , ff, f , /f} T - is vector, d= [d itj \, (i = 1, ...,4; j = 1, ...,4) - is the matrix 
of coefficients. The condition (17) represents the system of the linear equations on functions 
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bk(x), k — 1, 2, 3 and quantities d it j, (i = 1, ...,4; j = 1, ...,4). Solving this linear system of 
the equations we have: 

x 2 

b 3 (x) = d 13 ■ s ■ x - d 14 , ■ s—, (18) 

h(x) = d 13 ■ s (s + 1) + d 14 ■ x ■ s (l - |) , (19) 

h(x) = d 13 ■ (s) +du-(Y + s2 - s + l ) ( 2 °) 

s 2 

d n = du ■ (s 2 - s) + d u , d 2 i = d 13 ■ (s 2 + s) , d 22 = d u ■ — + d u , d 2 4 = Sd 13 , (21) 

3s 2 ( s 2 \ 

d-31 = d u ■ — , d 33 = -d 14 ■ — + s ) + d AA) d 4 i = 2d 13 ■ s 2 , (22) 



2 ' M \2 

s 2 

d A2 = d 14 ■ —, d 43 = d 13 ■ (s - s 2 ) , d 12 = d 23 = d 32 = d u = 0. (23) 

The existence of the non-trivial solution (18-20) means that the way by which extension 
(15) was made is suitable. 

Let us analyze the general form of the operator P 4 for the obtained solution: 

Pa = d 13 ■ s ■ Pi - ^|^P 4 2 + d u -d 13 -s + d 14 ■ (s 2 - s) (24) 

where 

^4 + < s + 1 >4^ w ;£ + < s - 2 >'4-*- < 25 > 

The operator P 4 (24) as well as operator P 2 (13), up to an additive constant , depends on 
two free parameters di 3 , d± 4 (24). From the general solution (24) we select two operators 
P4, P| (25) which do not depend on parameters d± 3 , du, rf 44 . We see that the coefficients 
of the operator P\ coincide with those of the operator P\ (14). This suggests the idea to 
implement the expression (25) for the operator P 2 1 (j V+1 ) for any N and subspace (15). The 
situation with the operator P^ is somewhat more complicated. The natural general guessed 
form for this operator reads 

d 2 d 

where the quantity depends on N and, besides, satisfies the conditions «o = 0, ot\ = 2 
as it follows from (14, 25). To find dependence a N from N, we shall consider the result of 
the action of the operator P 2 2 (7v+i) 011 the element /+ (15): 

AVi) (fn) =n(n + s-l-a N ).f++ ^ - ^ • f' +1 . (27) 
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If we put n = N in (27) and require that the subspace be finite-dimensional, the condition 
of cut off at n — N gives us a N = 2 • N, whence the operator -P 2 2 (at+i) looks like P^n+i) = 
x 2 -£z + (s — 2 • N) -x£ — x. As a result, we obtained a pair of operators that indeed leave the 
subspace M.2(n+i) invariant. Correspondingly, any linear combination P2(n+i) = Ki-P 2 \ 



2(JV+1) 

+ K2 • P 2 \n+i) wn ere k±, k 2 arbitrary constants has the same property. 

Following the way we described above, we found a series of concrete examples 1, 2, 4, 5, 
6. The list of obtained finite-dimensional subspaces together with the basic operators for 
which corresponding subspaces are invariant is given below. Meanwhile, the suggested way 
of the QES - extension (16) of the operators (4) is not unique. As the alternative approach, 
in an example 3 we used another "natural" procedure of QES - extension in that an integer 
n labeling the basis was added to the parameter of the function: a — > a + n. 

1) Finite dimensional function subspace TZ]^ = span{fQ, f£, f^} (N = 

0,1,2,....), dim (7^) = 2(N+1), formed by functions /+ = x n ■ F, [-;x] , f~ = x n ■ 
F 1 [s+^x] (n = 0, 1, .., N — 1, N) [16], is invariant for the operators = P 2 \n+i)i ^i = 



p2 

^2(JV+l)- 



Jr = *& + (* + i)£ 

J+ = x*£ + (s-2N).x£-x 



(28) 



dx 2 V ' dx 

The operators J±, Jf act on the functions fn( x ), fn(x) € Ti l N as follows: 



■n 



f+\ ( n(A n -l + s)-f+ + 2 -^. f~ 



+i 



fnJ \(n - s) (A n -l)-f~ + 8 (2B n - 1) • /+ 



j (f+\ = (f: + n(n + s). ft, + ^ • /" \ 
1 \fnJ \fn +n{n-s)- f~_, + 2ns ■ f+J 1 } 

, where A n = n — 2N, B n — n — N. The function f£ (x) satisfies the differential equation 
+ s £~ 1 ) fo (x) = 0. 
2) Finite dimensional function subspace 1Z 2 N = span{ft ■■■■> /jv> /cm ■■■■> /jv) (N = 
0,1,2, ....), dim(T^) = 2(N+1), formed by functions /+ = x n ■ ^ [ a g ;x] , f~ = x n ■ 
\P\ [s+i ; x ] ( n = 0' 1> ••' N —l,N) [16], is invariant for the operators J 2 , J 2 - 

J2^x£ + (1 + S-X)£ 

J 2 = ^£2 + (s - 2N - x) ■ x± + x (N - a) 
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The operators J% ', J 2 act on the functions /+(x), / n (x) G 7^ as follows: 

7+ f # ^ = T ^ - 1 + s) • / n + - 5n • + 2aB n /s ■ f~ +1 \ 
2 \fJ \(s-n)(l-A n ).f-+B n .f- +1 + s(2B n -l)-f.J 1 ] 

J = ( {a ~ n) f " + H {n + S) ' ti' 1 + a{1 + 2n) /S ' f A (32) 

2 \fnJ V (a + n + l)f-+n(n-s)-f-_ 1 + 2ns-f:_ 1 J 1 J 

, where A n = n — 2N, B n = n — N . The function (x) satisfies the differential equation 
(x£ + (s-x)£-a)f+(x) = 0. 

3) Finite dimensional function subspace 7Z% = span{fo, fi, /at} (N = 0,1,2,....), 
dim (11%) = N + 1, formed by functions f n = ^ [ a+ s n ;x] (n = 0, 1, .., N - 1, N) [16], is 
invariant for the operators J^: 



J+ = X 2JL + ( s ^ n - x) ■ x£ - ax 



(33) 

The operators J% act on the functions / n (x) G 1Z% as follows: 



Js (fn) = (n + a)-f n (34) 
J 3 (fn) = (sn +(a + n) C n ) ■ f n + (a + n) B n ■ f n+1 + n(a + n- s) ■ / n _i 

, where C n = N — 2n, B n = n — N . It is worthwhile to note that replacement 
(s — > s — 1, a — > a — | + y, N — > y — |) transforms the operators , J 2 into J^, 
J 3 ~ , if N is odd. 

4) Finite-dimensional function subspace 71% = span{fQ, f%, /q", f^} (N = 
0,1,2, ....), dim (7Z%) = 2(N+1), formed by functions /+ = x n ■ Airy (x) , /" = x n ■ 
-^Airy (x) (n = 0, 1, .., N — 1, N) [16], is invariant for the operators J^, J^: 



J- A =x^-(l + 2N)±-x^ 



' dx 2 yJ .^^yj dx 

J i = dx? - x 

The operators J^, act on the functions /^(x), fn( x ) e ^-n as follows: 



(35) 



J+ (fn\ = ( n{n-l)-f+_ 2 + 2n.f; i _ l \ 

" 4 \f-J {(1 + 2n) . /+ + n (n - 1) • f'J W 

T = ( {An ~ 2)n - #-1 + {2Bn ~ l) ■ f A (V) 

J *{fJ { 2B n .f: +1 + (A n -2)n.f-_ 1 J [6n 

where A n = n — 2N ', B n = n — N. The function /q 1 " (x) to satisfy the differential equation 
{& - x) fo + (x) = 0. 
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5) Finite dimensional function subspace 1Z% = span{fQ, f£, fo, • f N } (N = 
0,1,2,....), dim (7^) = 2 (N + 1), formed by functions /+ = x n ■ BesselK (u,x) , /~ = 
x n • BesselK (u + 1, x) (n — 0, 1, .., N — 1, N) [16], is invariant for the operators J$, J^: 

u 2 +u+x 2 



(38) 



The operators J^, J 5 act on the functions /+(x), f n (x) £ 71% as follows: 

7+ = ( ^ + n) (Z/ + ^ ' ^ + " ■ ^ ^ f oo^ 

5 \fJ \{n-l-v){A n -l-v).f--2B n .f:J 1 J 

7-W = ( ^( 1 + ^ + 2 ^)-/n- 1 -(l + 2n)-/- \ 
H/J V-(l + 2n)./+ + n(n-2i/-l)./-_J 1 J 

, where A n = n — 2AT, B n = n — N. The function /J" (x) satisfies the differential equation 

(^ + ^-(^ + ^ 2 ))/o + W=0. 

6) Finite dimensional function subspace 1Z% = span{fQ, f£, /^} (AT = 

0,1,2,....), dim (7^) = 2(N+1), formed by functions /+ = x n ■ ^ L/ 2 ;x 2 ] , /" = 
x n ■ 1 F 1 ; x 2 J (n = 0, 1, .., AT — 1, N) [16], is invariant for the operators Jq,J$: 

(41) 



6 dx 2 dx 



4 = (2^ 2 + l + 2N)£ + 2x(N- 2a) 

The operators Jg 1- , J 6 act on the functions /+(x), /~(x) G 72.^ as follows: 

6 \f-J V n(A n -2).f-_ 1 + (2B n -l).f+ + 2B n .f- +1 J 1 J 

7- (ti \ = ( 2 ( 2 « - ») • /n + " n) ' /»-2 + Son • /n-i \ ^ 
6 V/J V2n-/ n + _ 1 + 2(2o ; +l + n)./- + (n2-r i )./-_J 1 } 

where A n = n - 2AT, B n = n- N. 

For the proof of equalities (29, 31, 32, 34, 36, 37, 39, 40, 42, 43) it is enough to 
use rules of differentiation of special functions or their representations by the power series 
[16]. The commutation rules for the operators j£, J k (k = 1 . . . 6) are given in appendix 
A. Let us now consider an explicit example of application of the constructed subspaces to a 
physical system - two-photon Rabi Hamiltonian. 
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III. QES TWO-PHOTON RABI HAMILTONIAN 

The two-photon Rabi Hamiltonian (TPRH) is an obvious extension of the original Rabi 
Hamiltonian [18], which takes into account the atomic transitions induced by the absorption 
and emission of two photons rather than one [17]. The corresponding system as a whole 
is not integrable. Let us prove that (33, 34), at some choice of parameters, is invariant 
subspace for TPRH: 

H=^-a z + u.b + b + g (b 2 + (6+) 2 ) ■ (a + + a_) (44) 

where a z ,a± = a x ± i ■ a y are Pauli matrices, b and b + are the annihilation and creation 
operators respectively ([&, b + ] = 1). Following the approach of the article [17], apply the 
Bogoliubov transformation to the operators b and b + (b(t) = cos(t) ■ b + sin(t) ■ b + , b + (t) = 
cos(t) ■ b + — sin(t) ■ b ), where t G [0, 2ir] is a parameter. We will use the coherent state 
representation (Fock-Bargman representation): b = b + = z. The eignevalue problem for 
the operator (44) takes the form : 



(45) 



where a = g (b (t) 2 + b + (t) 2 ) = 2g(£ + z 2 ^, c = uu ■ b+(t)b(t) = % ■ 
- sin (2t) -£2 + 2 cos (2t) z-£ + sin (2t) z 2 + cos (2t) - lj . Removing the function ip 1 (z) 
from the system (45) we come to the fourth-order differential equation for the function 

ip 2 (z): 






L(i)V 2 (z) = 



1 f A A\ / A A\ 2E A CJn E 2 
(I -i ( I ( (I ''I—'''- - — 



ij 2 (z) = (46) 



uj \ /v / or e iU)" uj^ 

Let us try to reduce eq. (46) to the algebraic problem using the finite-dimensional pre- 
sentation of the operators J 3 (33). For this purpose, we shall present the operator L (t) 
(46) as polynomial on the operators J/, J% . To this end, we use a gauge transformation 
<p {z)~ x L (t) 4> ( z ) = P {J3 , J3) \x=(i{z)i where P (xi,x 2 ) is a polynomial of the second degree 
in non-commutative variables x\ and x^- We restrict ourselves by a special case of 4>{z) 
and fJ>(z): (ft (z) = (f) 1 (z) = exp (r/ ■ z 2 ) or <p (z) = 2 (z) = z ■ exp (77 • z 2 ), fJ>(z) = £ • z 2 . 
Thus, a search of a polynomial representation for L (t) is reduced to a search of coefficients 
of polynomial P (xi,X2) and constants 77, ^, t. The obtained solutions read: 
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Type I, s 

L(t )e" 
Type II, s = 3 

l(- /„: 

where t 



2' 4 2 ' 

2 



(2 (J 3 -) 2 + [J+ J 3 -] - 7J 3 + + 4J 3 - + d) 



(47) 



2' 4 2 " 



(2 (J3-) 2 + [j+, j 3 ] - 74 - 8J3- + c 2 ) u =52 



(48) 



^ circtcin f 23" 



Q w 



ii 



JV+1 
3\/2 



V2 /-< _ q "o nj2 N 1 1 (~< . 
— , - o^ja - - T + g, C 2 - 

N 2 + + f , £ = f = Two types of the solutions (47,48) have different 
symmetries with respect to the parity operator II = exp (mb + b) [17] ([II, H] =0). 

Further, we will consider the solution (47) (type-I) in more details, the second solu- 
tion (48) can be considered in a similar way Taking into account the equality (47), the 
problem (46) can be reduced to a search of vectors <p k satisfying the matrix equation 

( 2 {hY + [J3J3] ~ 7 33 + 4 J3~ + Ci) 0fc = °> where j£, h are finite-dimensional repre- 
sentations of the operators , J3" (33) in the subspace 1Z 3 N . For each of <p k) taking into 
account (47), we have L (t ) e v ' z2 <p k = that is equivalent to (46) under the condition 
■02 = ty\ — ^'^(Pk- The functions obtained are entire ones and belong to the Fock space: 
|V^)f = 0(|^exp(-^(^+^ 2 ))) (type I), |^)| 2 = O (\z\ p exp (z 2 + z 2 ))) 
(type II). We list the corresponding functions explicitly for the special case N = 2. 



type I: 



^2 



1^1 



1^1 

-1/4. 
1/2 ' 



-5/4. 
1/2 ' 



2a; 



5 



( 57 , 7yjA 
1^20 15 y 



+ 



_ _L V42 
3 ^ 30 



+ 1^1 



3/4. 
1/2' 




(49) 
(50) 

(51) 



type II: 



^2 



— rj-z 



1^1 



1^1 

5/4. > 2 

3/2' ^ Z 



1/4 " f 2 2 
31/2' ^ 



5 , Vio 
4 7 



(! + #)+ 1^1 



+ 

"9/ 
3/2 



(52) 
(53) 



12 



where £ = ^p, ?7 = § ? fj = V% — \- The second component ■0 1 of the spectral problem 
(45) can be obtained from the equation ip 1 = ^ • (^E + a — ■0 2 . In table 1 the values 
of relative frequencies f° r both types I and II are given for different dimensions of the 
subspace 1Z%. The obtained solutions (47, 48 and table 1) concern a non-resonant case 
^ 7^ 1 and are related to the known isolated solutions (Juddian solution), found earlier in 
[17]. It is worth stressing that in [17] the eigenf unctions were constructed on the basis of 
elementary function, whereas the solutions (47, 48) are constructed on the basis of functions 
1F1 [ a+ s n ; x] ■ For the values of the parameters s = 1/2, a = -1/4 - N/2 (type I), s = 3/2, 
a = 5/4 — N/2 (type II) the hypergeo metric function does not degenerate into polynomial 
(-a £ N). 

Tablel. Values of parameters Hamiltonian (44,45), for different d\m(7l%), ^ = 

E _ JV+1 _ 1 





3 


5 


6 


7 


8 


I 


2u 
u> 


0.44315 


1.68889 


3.03496 


2.72766 
3.60267 


2.10305 
3.74421 
3.90266 


II 


2u 


0.79838 


0.79838 


2.23006 
2.75234 


3.43545 


2.66128 
4.08801 


E 

U! 


1.23205 


2.38675 


2.96410 


3.54145 


4.11880 



IV. CONCLUSION 



To the best of my knowledge, the only previously known example of QES related to 
special functions was found in [19] .where these function appeared "by chance" for a partic- 
ular problem connected with quartic Bose Hamiltonians. Now, we developed a systematic 
approach of QES-extension that enables us to generate new QES operators based on spe- 
cial functions. This extends considerably the family of QES systems and can find physical 
applications, one of which (two-photon Rabi Hamilatonian) was discussed in the present ar- 
ticle. The main features of our approach include 1) the construction of the afline annihilator 
/C (.M2) [15]; 2) multiplication at the power function x n . One can think that the approach 
suggested in the given work, will give rise to further essential expansion of classes of quasi 
exactly solvable models. 
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VI. APPENDIX A. 

The constructed operators j£, (k = 1...6), satisfy to the following commutation 
relations: 

[j~,4]=S k (55) 

[4, s k ] = 4 ■ (Jk) 2 + 4 ■ 44 + 4 -s k + 4- j+ + 4 ■ r k + 4 (56) 

[J fc - S k ] = q • (J,") 2 + c 2 • (J+) 2 + c 5 " • J+ + c 6 • J fc " + (57) 

where A; = 1 ... 6, 4 — 1> • • • > are constants, theirs value are given in the table 2. The 
operators (55) have the general structure Sk = (3 3 (x) + (3 2 (x) + Pi (x) 4- + (3 (x), 
j3 m (x) G V n ijn = 0, 1, 2, 3). However, we do not list them here explicitly since the 
corresponding expressions are rather cumbersome. 

Table2. Values of constants (? = 1, . . . , 7) included in the commutation relations (55- 
57), C N = 2 + 2a + s, a N = 2 + N + 2a, (3 N = -4 • (1 + u 2 + 2N + u), 1n = (a-N) ■ 
(s + 1) • Cat, 5 n = -8 • (2a - N) ■ a N , S N = s ■ a ■ (s - N - 2), A N = s + N + 2a, 
B N = (2N-s)-(s-2- 2N), D N = (s + 1) • (s - 2 - 2N), G N = (a - N) ■ (s + 1). 
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fc = 1 


fc = 2 


k = 3 


k = 4 


fc = 5 


fc = 6 


ct 


o 


o 


o 


o 


o 


-6 


c7 
°i 


2 


2 


2 


o 


2 


o 




o 


o 


o 


6 


o 


o 




-4 


-4 


-4 


o 


-4 


o 




-2 


-2 


-2 


o 


-2 


o 


ct 


2 


( J AT 




-2 


o 


o 




o 


1 


1 


o 


4 


4 


4 


B N 


B N 


s-N 





1 -4iV 2 


12 + 32a 


C 6 


-2 


—Cn 


-A N 


2 












In 


Sn 








5 N 


C 7 





Gn 


s ■ a 





Pn 
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